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Abstract: 

The photonic and QCD radiative corrections at next-to-leading order are calculated 
for fermion-pair production at muon colliders, maintaining the full mass dependence and 
helicity information of the muons and the produced fermions. Higher-order effects of 
initial-state radiation are included at the leading logarithmic level. In the calculation 
particular attention is paid to the issue of gauge invariance in the vicinity of resonances. 
The most important corrections are presented in analytical form. The detailed numerical 
discussion concentrates on the corrections to the (s-channel) Higgs-boson resonances in 
the Standard Model and its minimal super symmetric extension. The results show that 
photonic initial- and QCD final-state corrections are very important in a precision study 
of Higgs resonances, but that (photonic) initial-final interferences are widely suppressed 
and only modify the non-resonant background. 
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1 Introduction 



The search for Higgs bosons will be among the most prominent ambitions in future 
high-energy collider experiments. Once one or even more Higgs bosons are found at the 
Tevatron [ |l| , the LHC [ , or a future e"'"e~ collider [ |[ , it is necessary to determine the 
properties of these particles, such as their masses, decay widths, and couplings, in order 
to reconstruct the scalar sector of the underlying spontaneously broken gauge theory. 
Current electroweak precision data [ Q and the direct search for a Higgs boson at LEP2 [ 
Q constrain the mass Mh of the Standard Model (SM) Higgs boson to Mh < 222 GeV and 
Mh > 114 GeV at the 95% C.L., respectively. In the Minimal Supersymmetric extensions 
of the SM (MSSM) the lower bound for the hght Higgs boson h is about Mh > 91 GeV [ 
H, while an upper bound of Mh < 130 GeV results from theoretical constraints (see e.g. 
Refs. [ 1^, 1^ and references therein). Thus, it seems very promising to search for a Higgs 
boson with mass ^ 140 GeV. Since the decay channel into on-shell W and Z bosons is not 
yet open in this mass range, such a Higgs boson has a very narrow width, rendering its 
direct measurement in a resonance distribution impossible at hadron and e"'"e~ colliders. 

A muon collider provides the unique opportunity to investigate Higgs-boson resonances 
in the s-channel within the clean environment of a lepton collider, since the Yukawa 
coupling of Higgs bosons to the initial-state muons is about 200 times larger than for e^, 
resulting in a relative enhancement by a factor 40000. Studies [ Q of the experimental 
feasibility of a muon collider show that the beam energy resolution can even be small 
enough to resolve a Higgs resonance with a width of some MeV, as it for instance appears 
for a SM-like Higgs boson with mass < 140 GeV. To be specific, for a SM Higgs boson 
of Mh ~ 110 GeV the relative accuracies achievable by scanning the /i"^/!^ — » H — >■ bb 
resonance with a luminosity of 0.2 fb^^ was estimated in Refs. [ |10|, ITT], |T2| to be ~ 1 — 3 x 



10~ and 0.2 for the mass and total width, respectively. Moreover, observing the reaction 
IJ,~^H~ — s> H ^ //, is complementary to Higgs production at e"'"e~ machines where the 
Higgs boson is produced via its couplings to Z or W bosons. The high accuracy at a muon 
collider also allows for a better discrimination between a SM and a SM-like Higgs boson in 
extended models like the MSSM as compared to hadron or e^e~ colliders. Furthermore, 
the existence of a SM-like Higgs boson h in the MSSM means that the two other MSSM 
Higgs bosons A and H are heavy; in this case a muon collider also proves to be the most 
promising machine to study the properties of these Higgs bosons. For instance, measuring 
the mass difference Mh — Ma yields an important check of the Higgs mass relations [ ^, 
in the MSSM, as pointed out in Ref. [ O. More information on the Higgs phenomenology 



at muon colliders can be found in the literature [ |T^, |TT|, |T2|, |T^, |T5|, |T6|, |T^ including 
in particular studies of CP violation [0, Higgs decays into r+r" pairs [|l5l, and flavour- 



changing neutral currents [|T^. In this context, the possibility to produce polarized muon 
beams proved to be very useful. 

In this paper we calculate the photonic and QCD corrections to the fermion-pair pro- 
duction processes /x'^/U^ —>■ ff in the SM and the MSSM. These corrections represent an 
important ingredient in a precision calculation of the Higgs resonances in these channels. 
In particular, photonic initial-state radiation (ISR) distorts the resonance shape, as it is, 
e.g., well known from studies of the Z resonance in e"'"e~ annihilation (see e.g. Ref. [ ^ 
and references therein). In contrast to the e^e~ case, however, we have to include the 
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finite-mass effects of botli initial-state and final-state fermions, since Higgs production 
proceeds via the so-called spin-0 channel which is mass suppressed. In detail, we calcu- 
late the full photonic 0{a) correction, which decomposes into ISR, final-state radiation 
(FSR) and initial-final interferences, and include the leading ISR effects beyond 0{a) in 
the structure-function approach [ |1^. In some existing studies ISR effects were already 
included at the leading logarithmic level as described in the appendix of Ref. [ 11|; how- 
ever, non-leading corrections have not yet been presented elsewhere. These effects are of 
particular importance if polarized muon beams are considered, since the leading logarith- 
mic approach ignores a possible spin flip of radiating muons, which mixes the spin-0 and 
spin-1 channels. QCD corrections are evaluated in next-to-leading order, i.e. in 0{as), 
also taking into account the full mass dependence. For b quarks in the final state the 
running MS mass is introduced. 

In our discussion of numerical results we mainly focus on bb final states and consider 
SM Higgs-boson resonances for Mh = 115 GeV and 150 GeV. In the MSSM we discuss 
two interesting scenarios, one where the H/A resonances lie on top of each other and one 
where they are sufficiently separated in mass. In the latter case, the masses are chosen 
large enough to enable an inclusion of tt production. Cross sections for polarized muons 
are only briefly discussed, although our calculation fully supports helicity eigenstate^ of 
the external fermions. 

The non-photonic, i.e. purely weak, electroweak radiative corrections of 0{G^) are not 
considered in this paper for the following reasons. It is known (see e.g. Ref. [|18[) that 
photonic corrections, in particular ISR effects, comprise a very important gauge-invariant 
class of radiative corrections to resonance processes of neutral particles. Another impor- 
tant source of corrections are renormalization effects, such as the running of the electro- 
magnetic coupling a. The latter universal effect is included by appropriately choosing 
the input value of a. Moreover, a fully convincing calculation of the non-photonic 0(0^) 
corrections is intrinsically connected to the Dyson summation of all relevant propagator 
corrections without losing gauge invariance. We discuss this issue in detail for the pho- 
tonic and QCD corrections and give an outlook to the non-photonic corrections, where the 
problem is more involved. In summary, our results can be viewed as a first step towards 
a precision calculation, and the weak corrections are part in one of the next steps. If 
one is only interested in a precise prediction of the Higgs resonance shapes, a reasonable 
approach would be to include the weak corrections to the Higgs Yukawa couplings on the 
Higgs resonances in the framework of effective couplings. This procedure maintains gauge 
invariance and the effective couplings can be easily combined with the photonic and QCD 
corrections discussed here. The most important weak corrections to the Higgs Yukawa 



couplings in the MSSM have already been discussed in Ref. [ 20 



The paper is organized as follows: In Section ^ we set our conventions and calculate 
the lowest-order cross sections. Section ^ contains our discussion of gauge invariance. The 
photonic corrections are treated in Section where our explicit analytical results for the 
ISR and FSR effects can be found. The QCD corrections are deduced from the photonic 
FSR results in Section |[ Section |] provides a detailed discussion of numerical results for 

^The generalization to other polarization states, such as transversely polarized beams, is straightfor- 
ward. 
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Figure 1: Lowest-order diagrams for /i ft^ f f with / 7^ /i 

the various corrections near the Higgs resonances. Our conclusions are given in Section 
and the Appendix provides some auxihary formulas. 

2 Notation and lowest-order cross sections 

We consider the process 

^^'{p,cr) + f^^ip',^') — ^ /(g,^) + /t^/^", (2.1) 

where the momenta p,. . .and helicities a,. . .of the muons and fermions / are given in 
parentheses. The heliticities take the values a = ±i, etc., but we will also take their 
sign to indicate the helicity. In lowest order the diagrams of Figure |l| contribute to the 
process ( |2.1|) . In order to treat the SM and MSSM in one go, we have abbreviated the 
Higgs bosons by where = H in the SM and = h, H,A in the MSSM. The field x 
denotes the would-be Goldstone partner of the Z boson; in the unitary gauge the diagram 
involving x exchange is absent. 

In the following we denote the vector and axial-vector couplings of the neutral gauge 
bosons to the fermions by 

^7/ = -<5/, a^f = 0, vzf = —Qf-^^^, az/ = -7/^, (2.2) 
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where I^j = ±i is the weak isospin of the left-handed part of the fermion field /, and 
Qf is its electric charge. The corresponding Feynman rule for the V^// vertex is given by 
ie7^(fv/ — cty/Ts), where e is the electric unit charge. In these couplings the parameters 
s„ and c„ are functions of the weak mixing angle 9^, 

sm 6'„ = s„ = 1 - c^, cv = -r^, (2.3) 

Ml 

and thus fixed by the ratio of the gauge-boson masses Mw and Mi. For the vertices 
5*// and Pjj of fermions with the neutral scalar and pseudo-scalar fields S and P we 
generically write the Feynman rules as \eg^f and egpfj^, respectively. For the pseudo- 
scalar would-be Goldstone field x the coupling factor reads 

9,f = -i^I^, (2.4) 

where m/ is the mass of the fermion /. Here we write Mz for the Z-boson mass instead 
of Mz for later convenience; the definition of Mz and its relation to Mz will be discussed 
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in Section || in detail. While the xff coupling is equal in the SM and the MSSM, the 
couplings of the fermions to the Higgs bosons, of course, differ in these models. We follow 
the conventions of Ref. 



21|1 also for the SM and define 



SM: 



MSSM: 



my 



9Hf 
9hu 

gnu 
9Au 



2c^s 



Mz' 



niy COS a 



Mz sin P 
1 TUu sin a 
2cvSw Mz sin (3 ' 

— cot p. 



2c„s„ Mz 



dhd 

9nd 

9Ad 



1 nid sin a 

Mz cos l3 
1 rrid cos a 
2cvSw Mz cos (3 ' 
I ma 

tan p. 



2cwSw Mz 



(2.5) 



where u and (i refer to up- and down-type fermions, respectively. In the SM the Yukawa 

coupling eguf is entirely determined by the ratio nif/v, where v = 1/ \J \plG^ is the 
vacuum expectation value of the Higgs field. In the MSSM the couplings additionally 
depend on the angles a and /3, which are related to the masses of the Higgs bosons. We 
follow the usual convention to use the mass Ma and tan/5 = f2/fi as independent input 
parameters, where f 1,2 are the vacuum expectation values of the two Higgs doublets in 
the MSSM. The angle a as well as the other Higgs masses Mh and Mh are then derived 
from Ma, tan/5, and the parameters of the gauge sector; to this end we employ the 
program JfDECAY [ ^ which calculates a, the dependent Higgs masses, and the Higgs 
decay widths F,^, thereby including also some leading radiative corrections [ to these 
relations. 

Sometimes it is more convenient to decompose the fermionic couplings into their right- 
and left-handed parts, which are proportional to the chirality projectors cj-i- = (1 ± 
75)72, instead of using vector, axial-vector, scalar, and pseudo-scalar couplings. Therefore, 
we introduce the right- and left-handed coupling factors for the various couplings 
introduced above, 



=F ay/, 



9% 



9s f, 



9pf 



(2.6) 



so that the Feynman rule for the V^^f f vertex is given by ie'y^{gyjU!+ + g^j^u^) and the 
rules for the scalar or pseudo-scalar field = S, P generically read ie((y'JjCo'+ + g~juj^). 
The Mandelstam variables for the process (|2.1|) are defined as usual. 



r\2 



t={p- qf 



u 



i\2 



(2.7) 



In the calculation of the tree- level amplitude A^o and of the one-loop amplitude M.i, 
which is described in the next section, we separate the fermion spinor chains by defining 
the following standard matrix elements (SME), 



(2.8) 
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with an obvious notation for the Dirac spinors Vfj_+, etc. The labels v and a in the 
SME refer to the insertions of 1 or 75 in the two spinor chains; for instance we have 
TVI™ = [v^+u^_][ufy5Vj]. The explicit expressions of the SME for all possible helicity 
configurations can be found in the Appendix. The SME represent 16 independent func- 
tions which form a basis for the tree-level and one-loop amplitudes, i.e. these amplitudes 
admit the representation 

■^n = E E F:y{s,t)Mr, n = 0,l, (2.9) 

i=l x,y=v,a 

with Lorentz invariant functions F!^\{s,t). The lowest-order amplitude reads 

Mo= -^v,o+ E -^P,o+ E -^3,0 (2.10) 

y=7,Z P=x,A S=h,H 



with 



e2 



S — My 

s - Ml 

Ms,o = -^^^MT, (2.11) 

where it is understood that the sums over P and S extend only over x and H, respectively, 
in the case of the SM. In a pure lowest-order prediction the propagator denominators 
(s — M^) lead to infinities in the resonance point s = M^. A perturbative description of 
these resonances requires a Dyson summation of propagator corrections that regularize 
the unphysical divergence. Here we signal the propagator modification, which is explained 
in the next section in detail, by putting carets over the masses in the propagators. For 
the photon we, of course, have = 0. 



Ignoring the Dyson summation for the moment, in Eq. ( |2.10|) the contribution of 
P = X results from the sum of x exchange and of the k^ki, part of Z-boson propagator in 
a general gauge, with k denoting the transferred momentum. The unitary gauge and 
the 't Hooft-Feynman gauge represent two interesting special cases. In the unitary gauge 
the X diagram is absent and the P = x part is completely due to the fc^/cy part of the Z 
propagator; in the 't Hooft-Feynman it is the other way round. 

Finally, the lowest-order cross section reads 

^o = ^^|^/dfiz E ^^(l + 2P_a)(l + 2P+a')|A^o(a,a',r,r')r, (2-12) 

where 



(7,a ,r,T 



/ Ami / 4m^ 

/?. = yi-^' ^/ = r-^' (2.13) 

are the velocities of the muons and produced fermions in the centre-of-mass (CM) frame, 
where also the solid angle f2/ of / is defined, and Ncj is the number of colours of /, which 
is 1 if / is a lepton and 3 if / is a quark. The degrees of polarization in the /i^ beams are 
denoted by P± and range from —1 to +1. 
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3 The issue of gauge invariance 



3.1 Preliminaries 

We first consider the case of a single resonance appearing in a propagator with mo- 
mentum transfer k and (renormalized) mass M. In lowest order the resonance appears 
for fc^ — > where the lowest-order propagator has a pole, and a proper description of 
the resonance requires a Dyson summation of the self-energy corrections to the resonant 
propagator. After this summation, in the vicinity of the resonance the amplitude Ai 
behaves as 

R 

^ ^ k^-M^ + S(A:2) + non-resonant terms, (3.1) 

where denotes the self-energy of the propagating particle. The resonance appears 

for a value of fc^ where k"^ — -|- Re{S(fc^)} = 0. Since S(fc^) develops a non-vanishing 
imaginary part in this resonance region, the propagator pole is shifted away from the real 
axis into the complex k^ plane. We denote the complex location of this pole by — IMF 
in the following. The real quantities M and T provide a consistent definition of the 
particle's mass and width, respectively. In particular, for gauge theories it was shown [ 
^ that this definition is related to the bare parameters of the theory in a gauge-invariant 
way. Making use of M and f , the naive propagator modification 

1 1 

(3.2) 



k'^-M^ P-M2 + iMr' 

which is known as fixed-width scheme, yields a possible description of the resonance in 
lowest order. 

Another definition of the particle's mass is provided by the zero of the real part of the 
resonance denominator in Eq. ( p.l| ), known as on-shell definition, which we identify with 
M in the following. This definition was, for instance, adopted for the Z-boson mass at 
LEPl. The quantities M and M coincide up to one-loop order, but differ at two loops 
and beyond. Making use of the usual renormalization conditions in the on-shell scheme, 
Re{S(M2)} = and Re{S'(M2)} = 0, the resonance behaviour of the propagator in 
Eq. ( p.l|) is obtained from the lowest order by the substitution 

11 . ^ 

(3.3) 



k2 _ M2 P - M2 + iMr(fc2) 

In the approximation of massless decay products the running width T{k'^) is related to 
the on-shell decay width F by r(fc2) = F x k'^ / M"^ 9{k'^). The parametrization ( |3.3| ) is 
known as running-width scheme. It can be checked easily that the parameter sets (M, F) 
and (M, F), when extracted upon fitting the parametrizations ( p. 21) and ( ^.3[ ) to data, are 
related by [ |2^ 



M = M/^l + 72 = M- — + ..., 7 = r/M, 

I F3 

r = F/Jl + 72 = F ^ + ..., 

/V ^ ' 2M2 

i? = i?/(l + i7) = /?(l-i^ + ...), (3.4) 



6 



where we have also included the relation of the corresponding residues R and R. 

At first sight it seems that the naive propagator modifications ( p.2|) and (p.3|) yield 
reasonable descriptions of the resonance at least in lowest order. However, this is not the 
case in general, since this modification potentially violates the gauge-invariance proper- 
ties of the full amplitude. Gauge invariance implies relations between self-energy, vertex, 
and box corrections (etc.) that in particular guarantee the cancellation of gauge- fixing 
parameters within the amplitude. These generalized Ward identities work order by or- 
der in perturbation theory, i.e. they are jeopardized by any kind of partial inclusion of 
higher-order corrections, such as the modification ( |3.2D which resums part of the self- 
energy corrections. Therefore, finite decay widths have to be introduced very carefully, in 
particular if radiative corrections are taken into account, which additionally complicate 
the situation. 

3.2 Survey of photonic and QCD corrections — gauge-invariant subsets of 
Feynman graphs 

Before we explicitly describe our treatment of the resonances occurring in the pro- 
cess it is useful to identify gauge- invariant subsets of Feynman graphs within the 
inspected class of corrections. 

At tree level (see Figure |ID only Z-boson and x exchange diagrams are linked by gauge 
invariance, while each diagram with photon or Higgs-boson exchange is gauge invariant 
in itself. This is also reflected by the fact that the Z- and Higgs-boson masses in the 
propagators can be viewed as independent parameters in the lowest-order diagrams. While 
this argument is trivially correct in the SM, where Mh is a free parameter, it works in 
the MSSM since the same diagrams appear, for instance, in a general two-Higgs-doublet 
model, where the masses are completely unrelated. Thus, the MSSM relations among the 
Higgs masses do not disturb this gauge-invariance property of the (tree-level) diagrams. 

Next we consider the photonic 0{a) corrections, which are induced by one- loop photon 
exchange diagrams and real emission of a single photon. This subclass of electroweak cor- 
rections is gauge-invariant, since a consistent U(l)eimg gauge theory of the neutral bosons 
(7, Z, and Higgs bosons) and all other fermions is conceivable, in which the process (|2.1|) 
proceeds in the same way as in the SM or MSSM in lowest order. If W bosons were 
involved at tree level this would be not the case. The photonic one-loop and bremsstrah- 
lung diagrams that are relevant in 0{a) are collected in Figures |^-^; the corresponding 
counterterm graphs are not shown explicitly. Of course, virtual and real corrections are 
separately gauge- invariant, since they correspond to different final states. However, these 
corrections are connected in the domains of soft and collinear photon emission, and these 
relations between one-loop diagrams and their real counterparts have to be preserved in 
the introduction of finite decay widths. Otherwise IR divergences and/or collinear singu- 
larities in the high-energy limit might not cancel properly. Therefore, we do not separate 
virtual and real corrections in the following discussion. 

Using the same argument as for the tree-level graphs, the diagrams of the photonic 
corrections can be further divided into gauge-invariant subsets according to the boson 
that is exchanged in the related tree-level graph. Thus, we have a gauge-invariant subset 
with internal photons only, one with Z-boson or x exchange, and one for each type of 
Higgs bosons. Furthermore, each of these subsets splits into three gauge-invariant parts 
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Figure 2: Diagrams for photonic vertex corrections to /i fi'^ ^ ff contributing to (a) 
initial-state and (b) final-state corrections with ip denoting any Higgs boson 
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Figure 3: Photonic box diagrams for /i fi^ ^ ff with V1V2 = 77, 'yZ, Z7 and ip denoting 
any Higgs boson 





etc. 



Figure 4: Lowest-order diagrams for /i fj,~^ ^ ff'y with ip denoting any Higgs boson. 
Graphs with the outgoing photon attached to the fermions /, / are not shown. 
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according to the charge factors of the two fermions to which the exchanged or emitted 
photon is attached. In this way we distinguish initial-state corrections, final-state correc- 
tions, and corrections connecting initial and final state. The relative correction factors are 
proportional to Q^, Q'j, and Qfj.Qf, respectively. The corresponding one-loop diagrams 
are shown in Figures 0(a), 0(b), and ^, respectively; the related squared diagrams of real 
photon emission should be clear. 

The 0{as) QCD corrections, which are only present if the fermion / is a quark, 
are closely related to the photonic final-state corrections of 0{a). The QCD diagrams 
are obtained from the photon exchange or emission diagrams upon replacing the loop- 
exchanged or emitted photon by a gluon. The QCD diagrams are, thus, classified exactly 
in the same way as the photonic final-state corrections. 

3.3 Treatment of Z-boson and Higgs resonances 

In the following we describe a simple way of introducing finite widths for the Z and 
Higgs bosons in the process ( |2.1| ) that retains gauge invariance in the presence of photonic 



and QCD corrections. 

As explained in the previous section, each Higgs-boson mass can be viewed as an 
independent parameter as long as we inspect only photonic and QCD corrections. Thus, 
a possible way of introducing finite Higgs-boson widths is provided by substituting the 
real Higgs masses squared by the complex pole positions in all Higgs propagators. Since 
all relations induced by gauge invariance, such as the cancellation of the gauge-parameter 
dependence in the photon propagator, are of algebraic nature, they are not destroyed 
by this substitution. For Higgs bosons the difference between barred and unbarred mass 
and width parameters will not be experimentally significant, so that there is no need 
to distinguish between barred and unbarred quantities any further. Thus, for the Higgs 
masses appearing in Eq. ( |2.1U| ) we set 

M^ = M2-iM^r^, (^ = h,H,A. (3.5) 

For Z-boson and x exchange the situation is much more involved, since the poles 
in the longitudinal part of the Z-boson propagator and the one in the x propagator are 
unphysical and have to cancel each other properly. We perform this cancellation explicitly 
in a general gauge, where the free propagators are given by 



V-^ ^ F - M| + iMzrz(F) l^^'^" F J P-ezMl 



with denoting the gauge-fixing parameter of the Z field. Here we have already taken 
into account that the transverse part of the Z propagator receives a contribution of the 
physical gauge-boson width after Dyson summation, while the longitudinal part and the x 
propagator do not. The actual form of the function Tzik"^) for the width is left open for the 
moment .0 In the diagrams for the tree level, for the one- loop, and for the bremsstrahlung 

^Contributions that do not change the resonance structure, such as real parts of self-energies or jZ 
mixing effects, are not included in the Dyson summation, since they are part of the purely weak correction 
not considered in this paper. 
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(a) k k (b) 7^ 7^ 



z 



Figure 5: Subdiagrams with Z and x exchange. The blobs represent tree-hke or loop 
parts, and the dots stand for external (on-shell) fermion or photon fields. 

corrections, shown in Figures ^]-|^, the propagators (|3.6|) appear in the two different ways 
schematically shown in Figure |^. The subdiagrams (a) contain Z or x fields on a tree-like 
line, i.e. cutting this line the diagram is decomposed into two disconnected parts; the 
subdiagrams (b) contain a Z or x line within a loop. As indicated in Figure ^ we denote 
the structures of the blobs in subdiagrams (a) by T'^'^ and [n = 1, 2). The functions T 
are nothing but (photonic parts of) amputated Green functions for off-shell Z or x fields 
contracted with the wave functions for the remaining external fields of /i"*", /, /, or 7 
which are all on shell. Since subdiagrams (a) form a subset of gauge-invariant diagrams 
(see previous section), the T functions obey the same well-known Ward identities (see e.g. 



Refs. [ |2^, |26[ ) as the full amputated Green functions that are on-shell contracted for all 



legs other than the Z or x line: 



- k^T^'^{-k) = -iM^T^i-k), k,T^'\k) = -iM^T^ik), (3.7) 

where ^k in the argument of Ti 2 are the incoming momenta of the subdiagrams. Using 
these identities the contribution Ai(a) of class (a) takes the form 

M^a) = T^''i-k)GfMT2^'''ik)+m-k)G^Hk)T,^{k) 

^ -iTf-^(-fc)T|^(fc) iT^i-k)T^{k) [k' + iM^T^jk')] 

F-M| + iMzrz(P) F[P-M| + iMzrz(F)] ' ^'^ 

where the arbitrary gauge parameter ^z has dropped out, as it should be. For diagrams of 
class (b) we proceed similarly. An important difference is, however, that the subdiagrams 
contained in the blobs do not obey Ward identities like Eq. ( |3.7| ) separately, because the 
blobs are linked. Instead the whole rest of the diagram without the Z or x line can be 
viewed as part of an amputated Green function with external on-shell particles /i"*", /, 
/ and two off-shell Z or x legs. These parts of Green functions, denoted by T^^'^'^(— /c, k) 
and T^^{—k, k), obey the Ward identity 

k^k.T^^'^-'i-k, k) = MiT^^{-k, k), (3.9) 

as we have also checked explicitly for the box diagrams. Maybe it is interesting to note 
that really the sum of all four relevant box diagrams on both sides of Eq. ( |3.9| ) is needed 
for the identity, i.e. it is not valid for the pairs of box diagrams obtained by interchanging 
Z and X lines. The contribution A^(6) of class (b) can, thus, be written as 

=/(^ {T^^'^''{-k,k)Gf^{k)+T^H-k,k)G^Hk)} 

■^Note that there is a global sign difference in the Ward identity in Ref. [ because of the different 
sign convention in the SU(2) gauge coupling e/sw 
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(27r)4 \A;2-M| + iMzrz(A;2) k"^ [k^ - + iMzTz{k^)] 

where the loop integration has been made exphcit. 

Now we have to specify how to choose the function Tzik"^) describing the Z-boson 
width. It is interesting to note that the simple choice of identifying Tzik"^) with the (on- 
shell) width Tz is not consistent, although Fz (M|) = Fz on resonance. The reason for this 
failure is an artificial infrared (IR) divergence in the loop integral of Ai{b), which appears in 
the second term in the curly bracket of Eq. ( p.lO|) for A;^ — > 0. The divergence is, however. 



avoided if the k'^ dependence of Fz(fc^) is taken into account, since Fz(A;^ — > 0) = 0. We, 
therefore, set Fz(fc^) = Fz x k'^/M^ which is exact for massless decay products and k"^ > 0. 
The case fc^ < 0, where Fz(/c^) actually vanishes, only occurs in the box loop integrals 
A4(b), and the effect that we do not set Fz(fc^) to zero there is of the order aFz/Mz, i.e. 
beyond our accuracy. The benefit of the chosen form for Fz(/c^) becomes obvious after 
inserting Fz(fc^) and reparametrizing Ai^a) and Ai^b) by the complex-pole parameters Mz 
and Fz: 



M(a) = (l + i7z) 



-1 



-iT^'^{~k)T,^Jk) + if^{-k)n{k) 
k'^-Ml + iMzFz 



where we have defined 7z = Fz/Mz and 



T^{k) = Jl + i7z T^{k), T^^-k, k) = {l + i7z) T^^{-k, k), n = 1,2. (3.12) 



The rescaling ( |3.12| ) is equivalent to substituting M| M| — iMzFz in and T^^, since 
these functions are proportional to M^^ and ^, respectively, because of the couphng 



of the X fisld to fermions. 
In summary, we have shown that the consistent use of the propagators ( |3.6| ) , where only 
the imaginary part, MzFz(/c^), of the transverse Z-boson self-energy is Dyson summed, 
together with the approximation Tzik'^) = Fz x fc^/Ml is equivalent to an overall re- 



placement Mz —>■ y — iMzFz in all graphs with Z-boson or x exchange before any 
Dyson summation is done, up to the global scaling factor (1 -|- i7z)~^- This observation 
again shows that no (algebraic) relations induced by gauge invariance are violated in this 
procedure. In our following calculation we make consistent use of the form ( 3.11| ), but 



drop the global scaling factor (1 + i7z)~^, since it influences squared amplitudes only at 
the level F|/M| or FzF<^/(MzM(p), which is beyond the accuracy of interest. This means 
that we set 

M| = m2 = M| - iMzfz (3.13) 

in Eqs. (|2j) and (|TT 



3.4 Outlook 

We emphasize that the procedure described in the previous section cannot be taken 
over without modifications to include the genuine weak 0{a) corrections (also denoted 
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by 0{G^) elsewhere) where the situation is more involved. The main complication is due 
to the fact that weak self-energy, vertex, and box corrections are all gauge-dependent 
and connected by gauge invariance. These relations have to be preserved by the Dyson 
summation. One possibility to introduce finite decay widths in the resonance propagators 
without losing gauge invariance is provided by the pole scheme [ ^ , where each propagator 
is expanded about its resonance pole and only the part of the residue, which is gauge- 
invariant, gets a finite width in the denominator. Details of the actual application of 
this scheme can also be found in Ref . [ ^ , where it was applied to the Higgs resonance 
in ZZ — > ZZ, including electroweak 0{a) corrections. In the pole expansions, of course, 
care has to be taken that the relations between photonic virtual and real corrections in 
the soft and collinear regions are still retained, which might additionally complicate the 
application. 

The discussion of the previous section raises the question about a possibility to consis- 
tently use complex mass values, which are derived from the complex locations — IMF 
of the propagator poles, as renormalized mass parameters. If the complex masses for the 
W and Z bosons are used in the definition ( |2.3| ) of the weak mixing angle, all generalized 
Ward identities, and thus gauge invariance, remains valid. This complex-mass scheme has 
been introduced in Ref. [ 3D| for lowest-order predictions. However, whether this scheme 
is also consistent for the treatment of radiative corrections is not yet clear, since the use 
of complex mass and coupling parameters is field-theoretically questionable. 

The field theoretically most desirable solution certainly consists in a consistent Dyson 
summation of all propagators (in fixed-order predictions), since this includes effects of 
running couplings directly in the propagator denominators. While this procedure in gen- 
eral violates gauge invariance in the conventionally quantized gauge theory, as explained 
also above, it was shown [ EG] that gauge invariance is preserved after Dyson summa- 



tion within the background-field formalism [ However, following this approach in 
a consistent one-loop calculation, the finite widths introduced via Dyson-summed one- 
loop self-energies are tree-level quantities, i.e. the prediction obtained this way is not yet 
fully (9(a)-corrected. Although it is not yet clear whether or how this problem can be 
solved consistently in a pragmatic way, the background-field approach seems to be very 
promising. 

The above considerations apply to a full inclusion of all weak 0{a) corrections both to 
resonant and non-resonant contributions. If one is only interested in a precise prediction of 
the Higgs resonance shapes, a simpler procedure is feasible. A reasonable approach would 
be to include the weak corrections to the Higgs Yukawa couplings in terms of effective 
couplings, i.e. the coupling parameters ( |2.5| ) would be dressed by the weak radiative 
corrections to the (on-shell) (pff vertices, and to treat the photonic and QCD corrections 
as described in this paper. This is possible for the following reasons. Firstly, it is known 
(see also the numerical results below) that the interference of the Higgs bosons with the 
continuous background is totally negligible, since the former results from the spin-0 and 
the latter from the spin-1 channel. Secondly, there is no problem with gauge invariance, 
because each set of Higgs-exchange diagrams is gauge-invariant at tree level and in the 
presence of photonic and QCD corrections, as explained above. 
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4 Photonic radiative corrections 



4.1 Virtual photonic corrections 



The virtual one-loop corrections add coherently to the lowest-order matrix element 
Aio given in Eq. (|2.10| ). Denoting the one-loop contribution to Aio that is induced by 
photon exchange by the one-loop-corrected squared matrix element reads 



\M\'' = \Mo\'' + 2Re{M^;'^M*o} + 



(4.1) 



where the dots stand for the genuine weak corrections of (9(6*^), not considered in this 
paper, and higher-order terms. The one-loop correction ai to the cross section is obtained 
by replacing |A^oP in Eq. (|2.12| ) by the one-loop contribution of As explained in 

the previous section, the one-loop correction M.i'^ receives contributions from vertex and 
box corrections, 

(4.2) 



(7) 



,(7) 

box) 



which will be considered in the next two subsections in detail. 

The actual calculation of the one-loop diagrams has been carried out using standard 

m 



and are evalu- 



techniques. The Feynman graphs have been generated with FeynArts 
ated in two independent ways. In the first calculation FeynCalc [ ^ is used to express the 
amplitudes in terms of SME and to reduce the occurring tensor integrals algebraically to 
scalar integrals with the Passarino-Veltman algorithm [ B3]. In the second calculation, also 



performed with Mathematica but without FeynCalc, we express the amplitudes in terms 
of SME and tensor coefficients. The reduction of the latter to scalar integrals is then done 
numerically. The results of the two calculations are in perfect numerical agreement (i.e. 
within machine accuracy for non-exceptional phase-space points). In both calculations the 
scalar integrals are evaluated using the methods and results of Refs. [ |3j, where 



ultraviolet divergences are regulated dimensionally and IR divergences with infinitesimal 
photon and gluon masses, and mg, respectively. The renormalization is carried out in 
the on-shell renormalization scheme, as e.g. described in Ref. [ p6l . 



4.1.1 Vertex corrections 

In the on-shell renormalization scheme the external self-energy corrections are ex- 
actly cancelled by counterterms. Thus, the virtual vertex corrections induced by photon 
exchange are entirely given by the diagrams shown in Figure ^ and the counterterm con- 
tribution for the vertices. Among the relevant renormalization constants, only the one for 
the masses, and 6m f, and for the fields, 6Z^ and SZ^^ , receive photonic corrections, 
which can be easily derived from the general results given in Ref. [ ^ . We immediately 
give results for the renormalized vertices. 

The vertex functions for the fermionic couplings of a general vector boson V, of a 
pseudo-scalar boson P, and of a scalar boson S are parametrized by 

^yiQ^Q) = ie7p(w/ - av/75) 

7p(w/ - av/75)Ai(s,m/) + av/7p75A2(s, m/) 
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2mf{q-q')p 2mf{q + q')p 

^A3(s,m/) + ayf ^75A4(s,m/) 



An 



T^^{q,q') = iegsf h + Q) ^K^{s , m j] 



(4.3) 



where the momenta q and q' of the vertex functions F correspond to the outgoing fermion 
and antifermion, respectively, and s = {q + q'Y- Note that we have aheady assumed that 
the fermion fields are on shell, i.e. we have used q^ = q'"^ = m'j and terms that vanish 
after contraction with Dirac spinors have been omitted. The virtual corrections to these 
vertices are parametrized by the functions Aj(s,m). At the one-loop level these are given 
by 



Ai(s, m 
A2(s,m 

A4(s, m 
A5(s,m 

Afi('s,m) 



-3 [Bo{s, m, m) — Bq(0, m, m) — 2] — 2(s — 2m^)Co(m^, s, m^, m^, m, m) 
f Am^B'^iw? , m-y, m). 



5m 



s — Am? 
s 



s — Am? 

3s - An? 



[i?o(s, m, m) — -Bo(0, m, m) — 2] , 
[i?o(s, m, m) — -Bo(0, m) — 2] , 

[-Bo(s, m, m) — -Bo(0, m, m) — 2] 



s — Am? 
Ai(s, m) — A2(s, m) — K^i^s, m) 

2 — 2(s — 2m^)Co(m^, s, m^, m^, m, m) + Am^B'^iw?^ m^, m), 
2 : — - [Bq{s, m, m) — Bq{0, m, m) — 2\ 



s — Am^ 

— 2(s — 2m^)Co(m^, s, m^, m^, m, m) + Am'^B'^^m'^, m^, m). 
The appearing scalar integrals are easily calculated to 
5o(s, m, m) - fio(0, m, m) - 2 = ln(xj. 



(4.4) 



Co(m'^ , s,m'^ ,m^,m,m) = — — 



ln(^)ln(x,) + ^ln2(x. 



Km^ 

^2 



2Li2(l + x, 



TT 



27riln(l + Xs 



Bo{m'',m^,m) 



dBo, 2 ^ 



1 




hlnf^^Y 




1 - 






\m J . 



where 



/ Am^ 
1 + ie. 



X a 



/9s -1 



(4.5) 



(4.6) 



and Li2(x) denotes the usual dilogarithm. As indicated in Eq. ( ^4.4| ), not all of the functions 
Aj are independent. Owing to the Ward identity (|3.71 ), the correction to the xff vertex, 
and thus the function A5, is completely fixed by the correction to the Z// vertex. 
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When applied to the process ( p.l|) , the parametrization (|4.3| ) can be directly inserted 
to calculate the final-state corrections, while the corrections to the initial-state vertices 
follow from obvious substitutions. The complete vertex correction to the amplitude is 
given by 



vprt 



a 



\ [QlAi{s,m^) + Q}A,{s,mf) 



V=7,Z 
+ 



M 



v,o 



V 



Qlay^A2{s,m^)(vyfM1'' - ayjMI 
+ Q}aYfA2is,mf){vY^M'i'' - ay^MI 



+ 



a 
At: 



Anii 
s 



QlA^{s,m^) + QjA^{s,mf) Al^ 



J2 [QlA5{s,m^) + Q}A5{s,mf) 



P=X,A 



+ ^ E [QlMs,m^) + Q}Ms,mf) 



Ms,o, 



(4.7) 



S=h,H 



where we made use of the lowest-order structures A^v,o, -Mpfi, and A^s,o defined in 
Eq. (|TTD. 



4- 1-2 Box corrections 

Compared to the vertex corrections, the calculation of the box diagrams, which are 
shown in Figure ^ involves quite a lot of algebra. In fact, the full result for the box 
corrections is too lengthy and involved to be reported here. Instead we give the most 
important gauge- invariant part of the boxes, denoted by Albolres) which contains the 
complete IR structure and all resonant contributions. To this end, we split the box 
correction according to 

^Sx = -M£U + '^Sx,conf (4.8) 

The finite remainder A^bolcont; which describes box corrections to the continuum, will 
not be given analytically, but it is included in the numerical results below. 
The part A^S.res 

is obtained from the full box diagrams upon setting the integra- 
tion momentum in the numerator of the diagram to the value that corresponds to zero- 
momentum transfer by the photon. If two photons are exchanged in the box, there are 
two such values of the integration momentum, leading to two independent contributions. 
Obviously, this procedure extracts the complete IR-sensitive part of the box diagrams, as 
can be shown by simple power counting: if is the photon momentum, the difference 
between the full box and the extracted part contains a factor q^, which cancels the (loga- 
rithmic) IR divergence (g^ 0) of the denominator. The integration domain near q^ —* 
is, however, also the only source of potentially resonant contributions of the box diagrams. 
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A resonance requires s ~ M]^, where X denotes any boson exchanged at tree level, and 
the only way to preserve the condition ~ s ~ M]^ {qx denoting the momentum of 
X) in the box diagram is to require that the photon, which also connects initial and final 
state, is soft {q^ —>■ 0). Otherwise the X propagator becomes non resonant. Since we did 
not set — i> in the X propagator denominator in our definition of box res 5 ^^e full 
leading resonance behaviour of each box is still contained in this part of the box ampli- 
tudes. The actual calculation of A^bolres rather simple, since it involves only scalar 
integrals, and the result is 



•^bol.res - Qt^Qf- 
TX 



E 



[ml + m}- t)Do{t, Mx) - {t ^ u) 



X=7,Z,h,H,x,A 



X {s-Ml)Mxfl/{l + 5^ 



(4.9) 



where Sx^ = 1 for X = 7 and Sx^y = otherwise. The gauge invariance of this part of 
the one-loop amplitude simply follows from the gauge invariance of the Q^Q / part of the 
IR divergence, which is completely contained in A^boLres- The IR-divergent scalar box 
integrals Dq, including the generalization to complex masses Mx, have been calculated 
in Ref. f|35|l; the results read 



Do{t,Mx) = Do{ml,ml,mj,m'j,s,t,my,m^,Mx,mf) 



M'x] 



A(t,m2,m^) 



2\n{xt) 



ln(l - xl 



In 



p,(T=±l 



^ [Li2(xtx^x^) + ( \Yi{xt) + ln(xp + In(x^)) ln(l - XtX^^xf, 

for X ^ 7, 



+ y + U2{xl) + ln=^(x^) + \n^{xf 



Do{ml, ml, mj, mj, s, t, m^, m^, m^,mf) 

1 / \ 

21n(a;t)ln I — ^ + 



s^X{t,ml,m}) 
where we have used the usual definition 

X{x, y, z) = x^ + + — 2xy — 2xz — 2yz 

and the abbreviations 

A-1 



(4.10) 



(4.11) 



Xf] 



Pt + r 



'1 



Am^mf 
t - [m^ - mf)' 



+ ie, 



aX 



1 + A 



aX 



\ 



1 



4m^ 



(4.12) 



Note that the explicit factor (s — in Eq. ( [1.9|) is cancelled by the pref actor in the 
functions, but the resonance factor in the lowest-order structure M.xfi survives. 
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It is instructive to inspect the part A^boxrcs from a different point of view. Eacli 
contribution from X 7^ 7 can also be viewed as tlie leading term in an expansion of the 
corresponding box diagram about the resonance pole. Hence, these terms are nothing 
but the virtual part of the so-called non-factorizable corrections which were discussed in 
detail in Ref. [ ^ for several processes with W and Z resonances.^ From the results given 
there it is also clear that the resonant part of A^boires ^^^^ ^e completely compensated by 
its counterpart from soft photon emission. Thus, we can expect that the sum of virtual 
and real initial-final corrections to the resonances will be very small, as it is also known 
(see e.g. Ref. [ |18| and references therein) for the Z resonance in e"'"e~ collisions for a long 
time. 

Since the continuum contribution A^boLcont is gauge- invariant and does not involve 
resonant terms, there is no need to introduce complex masses in this part. Thus, we keep 
all masses real in the tensor integrals there. Introducing complex masses would change 
the cross section by terms of 0{arx/Mx) which is beyond the desired level of accuracy. 

If one is only interested in resonance cross sections, a simple approximation for the 
corrections in the vicinity of the resonance at s ~ Mx is provided by the neglect of the 
continuum contribution A^bol cont- This approach is known as pole approximation and, for 
instance, successfully applied for Drell-Yan-like W production at hadron colliders [ E3] 



and W-pair production in e"'"e annihilation [ ^ . One advantage of this approximation 



is certainly its simplicity compared to the full off-shell calculation. 



4.2 Real photonic corrections 

4-2.1 Amplitudes for fi~ fi^ — > ff'j 

The real photonic 0{a) corrections are induced by the process 

fi~ip,a)+fx+ip',a') /(g,r) + /(g',r')+7(^,A) (4.13) 

in lowest order, where k and A denote the momenta and helicity of the emitted photon, 
respectively. Typical diagrams for this process are shown in Figure ^ For the case of 
the SM the helicity amplitudes have already been given in Ref. [ 0. In the following we 
modify and generalize these amplitudes for our purposes. For completeness we repeat here 
the basic definitions needed to evaluate the amplitudes, but suppress most of the details 
concerning their derivation, which was done within the Weyl-van der Waerden (WvdW) 
spinor technique. We consistently use the variant of the WvdW formalism described in 
Ref. [ 0], where also references to other spinor techniques are given. 

In the WvdW formalism, each of the momenta p'^, etc., is written as a complex 2x2 
matrix etc., consistently denoted by the respective capital letter in the following. 
These matrices are decomposed into Weyl spinors as follows, 

Pab=T. Q%=Y. f^Af^B. K^B = k^kB, (4.14) 

i=l,2 ' i=l,2 

^Actually M^^J^ ^.^^ differs from the non-factorizable corrections by non-rcsonant contributions, since 
the off-shelhiess s is set to its resonance value in the latter. 
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where a dotted spinor index indicates complex conjugation. Using polar coordinates for 
the momentum p, 



= {po, |p| cos0psin6'p, |p| sin0psin^p, |p| cos6'p), 
the (covariant) spinors Ki^A are given by 



Po + iPl 



e-'*f cos f 
sin % 



«2M 



Po - iPl 



sm 



COS ■ 



Moreover, contravariant spinors are defined according to 



+1 
-1 



(4.15) 



(4.16) 



(4.17) 



The application to k[ ^ and pf]^ is obvious. The spinor /c^ is constructed from the mo- 
mentum in the same way as Ki^a from (the counterpart of K2,a vanishes for light-like 
momenta). For two Weyl spinors and ipA the antisymmetric product 



= 0aV^^ = M2 - 02^^!, (^V')* = = (01^2 - 02V^l)' 



(4.18) 



turns out to be very useful, since both Minkowski products and Dirac spinor chains can be 
written in terms of such products. For helicity eigenstates of fermions, the Dirac spinors 
can be expressed in terms of the associated momentum spinors in a very simple way. In 
our case we can write the Dirac spinors as 



u, 



A ' 



I A I ' 



(4.19) 



with the actual insertions 

(0,V^) = 





-/«2) 


for a = 


+, 


(/«2 


Ki) 


for 0" = 






-P2) 


for r = 


+ , 


(P2, 


Pi) 


for T = 


5 



(0', V^') 





for 


a' 


( — K2, K.[) 


for 


a' 




for 


t' 


(-p'2,P'l) 


for 


t' 



(4.20) 



for the various polarizations. 

The helicity amplitudes M.^ for the process ( |4.1^ ) read 



;W^(cT,(T',r,r',A) = V2e^ 



(4.21) 



LV=7,Z (p=x,h,H,A 

where the two generic functions A^^'^ and A^^^ describe spin-1 and spin-0 exchange, re- 
spectively. In the SM the sum over Lp extends only over H and x- For A = +1 these 
functions are given by 



4^)(a,a',r,r') 



Q^ikP'Pk) 



QfikQ'Qk) 



2{p ■ k){p' ■ k)[{q + q'Y - M2] 2(g ■ k){q' ■ k)[{p + p'f - M^] 
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{q + q'f - M\ 



V 



(p' ■ k) 



{9^fm*m+9yf{kiY{H 



Qf 



V 



9Yf{kvT 
{q' ■ k) 



Qf{kQ'Qk) 



{9Uk<pr{<Pv)+9y,{k^pr{^P'v) 

Q^jkP'Pk) 

4(g ■ k){q' ■ k)[{p + p'y - M2] 4(p ■ A;)(p' ■ + q'^ - M^] 

X + 9^M'^r) {9;f{^'v) + 9;f{^vr) 

Qu9r„„{k(p')*{ki^)* / ^ , , \ r 1 i 



Qf9-f{k0*{kvr 



9UH') + 9;;M'i^y 



2{p ■ k) 2{p' ■ k) 

1 1 

+ 



2{q-k) 2{q'-k) 



(4.22) 



Here we made consistent use of the chiral couplings introduced in Eq. (|2.6|) , and the dots 
stand for ordinary products of four-vectors. Moreover, the following abbreviations are 
used, 

{kp'pk)= y: {kK',r{K,K[){K,k)\ {kQ'Qk)= y: {kp'd*{pjp[){pjkr. (4.23) 

The amplitudes for A = — 1 follow from discrete symmetries. For instance, a parity 
transformation leads to the relations 



{9±)*^9j.,M*.^M.. 



(4.24) 



M^{-cr, -a', -T, -t', -A) = sgn{acx'TT')M^{cr, a', r, r'. A)' 
and CP implies 

= -sgn(^^Vr')^,(a, a', r, r', Xy 

(4.25) 

The substitutions for (g^)* and M* ensure (in addition to the interchange of and g 
in the P transformation) that complex couplings and masses are not subject of taking the 
complex conjugate of the amplitude. 

The bremsstrahlung contribution cr^ to the cross section is obtained from the transition 
amplitudes Ai-y by the phase-space integration 



^/d<^'7 E \il + 2P^a)il + 2PW)\M,ia,a',r,r',X)\' 



(4.26) 
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where the phase-space integral is defined by 



/ 



d$. 



f d^q 


f d^q' f 


/ (27r)32go^ 


(2vr)32g^ J 



d^k 



(27r)32A;o 



{2Try 5{p + p -q-q -k). (4.27) 



^.^.^ Treatment of infrared singularities 



The phase-space integral ( [4.26|) of the real photonic corrections diverges in the IR 
region, i.e. if the photon energy feg tends to zero. As done for the virtual corrections, we 
regularize this IR divergence by an infinitesimal photon mass m^, i.e. the integral (|426|) 
has to be modified appropriately. We have applied two independent methods for this task, 
known as phase-space slicing and subtraction. 

In the phase-space slicing approach, the soft-photon region < ko < AE is excluded 
from the integral (|4.26|) by the energy cut AE. If AE is chosen sufficiently small the 
photon-emission cross section can be treated in the soft-photon approximation, where 
it factorizes into the differential lowest-order cross section duo and a universal eikonal 
factor (see, e.g., Ref. [ ^). Integration over k in the soft-photon region yields the simple 
correction factor (5soft to dao, 



4oft = - {QlfiP,) + Q'ffiPf) + Q,Qf m - g{u)]} 

where the auxiliary functions / and g read 



(4.28) 



9it) 




2{ml + m) -t) 



^X{t,ml,m} 



(4.29) 



with the abbreviations 



at 



^\{t,ml,m 



2ml 



Vt 



2{alml-m)) 
siat - 1) 



(4.30) 



and the analogous quantities a„ and Vu-, obtained from at and f t by t ^ u. It can be 
checked easily that all IR-singular logarithms Inm^ of 5soft are compensated by their 
counterparts of the virtual corrections. When adding the parts of cr^ obtained from the 
soft and hard regions, the logarithms of the auxiliary cut AE cancel in the limit AE —* 0. 
Note that in the hard regime these logarithms result from a numerical integration. The 
cancellation proceeds up to terms of 0{AE/e), where e is a typical scale in the process. 
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Since e ~ is of the order of some MeV in the vicinity of the hght Higgs resonances, the 
cut AE has to be chosen extremely small there in order to achieve reasonable accuracy 
in the phase-space integration of cr^. 

This numerical subtlety is circumvented in subtraction methods. Here, the general 
idea is to subtract and to add a simple auxiliary function from the singular integrand. 
This auxiliary function has to be chosen such that it cancels all singularities of the original 
integrand so that the phase-space integration of the difference can be performed numeri- 
cally. Moreover, the auxiliary function has to be simple enough so that it can be integrated 
over the singular regions analytically, when the subtracted contribution is added again. 

We have applied the subtraction method presented in Ref . [ ET| , where the so-called 



dipole formalism, originally introduced by Catani and Seymour within massless QCD, 
was applied to photon radiation and generalized to massive fermions. Since the actual 
implementation of this formalism in our calculation does not involve any subtleties, we 
do not describe it here but refer to Ref. [ ^ for details. In particular, the application to 



the process — > z^e^^e, which receives the same initial-state corrections as the general 
// production process (^.l]), is spelled out there explicitly. 

The numerical results shown below have been obtained using the subtraction approach, 
since the corresponding Monte Carlo integration errors are considerably smaller than for 
the slicing method if the same number of integration points is used in either case. We 
have checked that the results of the two methods agree within the integration errors. 



4-2.3 Leading initial-state radiation beyond 0{a) 

The emission of photons coUinear to incoming high-energetic charged fermions, such 
as e^ or /i^ with energies of at least several 10 GeV, leads to corrections that are en- 
hanced by large logarithms. These leading logarithms are process independent, and their 
contribution to the cross section can be obtained by a convolution of the lowest-order 
cross section with a structure function, which is explicitly known up to 0{a^) [ |19|, ^ . 



Following the conventions of Ref. [ 43 1 , this convolution reads 



dx (^(2a, X, Q^) <7q{xs) 



(4.31) 



where o"o(xs) is the lowest-order cross section at the reduced CM energy squared xs. 
The flux function (f){2a, x, Q^) accounts for radiation from both initial-state muons and is 
connected to the QED splitting function F^^, which is up to order 0{a^) given by 



x,Q'^) = (j){a,x,Q'^) 

_ exp (-|Asr7e + |Asr) /?] 



r (l + |/3isr) 
PIk r 1 + 3x2 



1 — X) 2 



/5] 



ISR , 



l + x) 



aSR 

32 

ISR 



( 1 + 6x 1 
I— ln(x) + 4(1 + x) ln(l - x) + 5 + x j 



(3' 



384 



'l + x) 



6Li2(x) + 121n^(l 



X 



37r^ 



;i + 8x + 3x2) ^ g(^^ + 5)(1 - x) ln(l - x) 
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+ ^(39 - 24x - Ibx^) 
where /3isr contains the leading logarithm (LL), 



+ 12(1 + x^) ln(x) ln(l - x) - ^(1 + 7x^) \yi^{x) 
1 



(4.32) 



/3: 



ISR 



IT 



(4.33) 



Note that the scale is not fixed within LL approximation, but has to be set to a 
typical scale of the underlying process; for the numerics we use = s. In Eq. ( [4.32| ) 7e 
is the Euler constant and T{y) the gamma function, which should not be confused with 
the splitting functions. Note that some non-leading terms are incorporated, taking into 
account the fact that the residue of the soft-photon pole is proportional to ln(. . .) — 1 
rather than the logarithm for the initial-state photon radiation. 

We add the cross section ( |4.31| ) to the one-loop result and subtract the lowest-order 
and one-loop contributions ctll,! already contained within this formula. 



0"LL,1 



dx 



(5(1 -x) + 2r^^'^(x,g') ao{xs), 



(4.34) 



in order to avoid double counting. The one-loop contribution to the structure function 
reads 



r 



LL,1, 



x,Q' 



/3, 



ISR 



l+X^ 



X 



AsR ,. 

iim 

4 e^O 



5(1 - x) ( ^ + 2 Ine) + ^^(1 - X - e)\^ 
2 / 1 — X 



In summary, the radiatively corrected cross section a reads 



a = ao + a'f^ + + ((Tll - cxll,i) 



(4.35) 



(4.36) 



where al^'^ is the virtual photonic one-loop correction. Note that the uncertainty that is 
connected with the choice of enters now in (9(a^), since all photonic 0{a) corrections, 
including constant terms, are taken into account. 



5 QCD radiative corrections 



As already explained in Section the diagrams for the 0{as) QCD corrections 
are obtained from the diagrams for the photonic final-state corrections upon replacing 
the loop-exchanged or emitted photon by a gluon, provided / is a quark. Simple colour 
algebra shows that the photonic final-state corrections turn into the QCD corrections upon 
replacing the coupling factor Q'ja by CpUs, where Cp = 4/3 is the quadratic Casimir 
operator of the quark representation. More precisely, this substitution applies to the 

to the squared matrix element 



matrix element AifJrt for the vertex correction in Eq. 



|A^^P for hard-photon emission obtained from Eq. ( |4.21| ), and to the soft-photon factor 
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(5soft of Eq. ( ^4.28|) . The contributions proportional to or Q^Qf, of course, have to be 
set to zero there. We always evaluate the running strong coupling constant a^ifi) at the 
two-loop level including five active flavours in the running, i.e. the top quark is decoupled 
from the running of as(/i). The renormalization scale is set to fi = ^/s, which is the total 
CM energy. 

The case that / is a b quark is of particular importance phenomenologically. For this 
case it is known that the Higgs Yukawa couplings iphh receive very large QCD corrections 
of the form asln(mb/MH) that are connected with the on-shell mass renormalization of 
the b quarkQ. The convergence of the perturbative series is improved considerably after 
introducing the MS running b mass mb(/i) at an appropriate scale /x. To this end, we 
rescale the iphh coupling g^]^ in all matrix elements according to 



mb 



h,H,A, 



(5.1) 



which replaces the on-shell mass mb by the running b mass mb {\fs) renormalized at the 
scale /i = v^. Note that /i^ = s corresponds to the virtuality of the s-channel Higgs 
boson ip. The running b mass is evaluated following the results of Ref. [ ^ (which are 
also used in HDECAY) . In order to avoid double counting, the 0{as) contribution of the 
replacement ( p.l| ) has to be subtracted from the one- loop matrix element Aii. Using 



mb(Ai) = mb |l - 
the consistent modification reads 

(a/s) 



as(yu) 



Ml Ml 



+ 



TT 



TT 



In 



\mi 3 



+ . . . 



(5.2) 



+ 



mt 



J ip=h,R,A 



g^b^g^bis) 



(5.3) 



Finally, we remark that we have compared our results against the known analytic 



form of the inclusive next-to-leading order QCD corrections [ ^ to the on-shell decays of 
scalar and pseudo-scalar Higgs bosons into heavy quarks. Switching off all non-resonant 
diagrams in the calculation of the QCD corrections to the resonance cross sections, the 
relative QCD corrections to the peak cross sections reduce to the corrections to the cor- 
responding partial Higgs decay widths. We find agreement between our numerical results 
for the relative corrections to the peak cross sections and the corresponding analytical 
results for the partial Higgs decays within the Monte Carlo integration errors. 



6 Numerical results 



^Actually, analogous logarithms are present in the photonic FSR corrections, but there is no need for 
a special treatment of those corrections, since they scale with the much smaller factor Q^a instead of 
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6.1 Input parameters 



For the numerical evaluation we use the following set of SM parameters [ ^| 



Mw 
Mz 

m,(lGeV) 



1.16639 X lO^^GeV" 
0.1181, 
80.419 GeV, 
91.1882 GeV, 
105.658357 Me V, 
4.62 GeV, 
190 MeV, 



a(0) = 1/137.0359895, 



171^- 

rrir. 



2.12 GeV, 
2.4952 GeV, 
1.77703 GeV, 
174.3 GeV, 
1.42 GeV. 



(6.1) 



Note that part of the input ( |6.1| ) is only needed for the calculation of Higgs-boson param- 
eters for which HDECAY [ ^ is used; in particular, this is the case for ms(l GeV) which 
is the MS running strange mass at the renormalization scale 1 GeV. 

In order to absorb part of the renormalization effects in the electroweak couplings, 
such as the running of a((5^) from = to a high-energy scale and some universal 
effects of the p parameter, we derive the electromagnetic coupling a = e^/(47r) from the 
Fermi constant according to 
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(6.2) 



when calculating the couplings for the lowest-order cross sections. In the relative photonic 
corrections we use a(0) as coupling parameter, which is the correct effective coupling for 
real photon emission. Thus, the 0(a)-corrected cross section scales like a{0)aQ^. In the 
following both lowest-order and corrected cross sections for bb are consistently 

evaluated with a running b-quark mass, as described in Section |[ 

In the MSSM, we take tan (3 and the mass Ma as the basic input parameters of the 
Higgs sector. For the evaluation of the Higgs decay widths we need additional input 
parameters, which are chosen as 



/i = 300 GeV, M2 = 200 GeV, 



M 



II 



Mr, 



At = Ai, = 1.5 TeV, 



M-, 



M 



Ll 



Mr 



Mr = Mt = 1 TeV, 



(6.3) 



where the notation follows the description of HDECAY [ 22 



Table |I] contains the results, as obtained from HDECAY^ for the Higgs-boson masses 
and decay widths, which are used in the following SM and MSSM evaluations. 



6.2 The process /x 



-,,+ 



bb in the SM 



We start our numerical discussion by considering the polarized cross sections for 
— > bb in the SM for Mh = 115 GeV. The l.h.s. of Figure ^ shows the lowest-order 
cross sections and their (photonic and QCD) corrected counterparts. We first concen- 
trate on the lowest-order cross sections, where the Higgs-boson resonance occurs only in 
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SM 


MnlGeV] 


115 


150 




rnlGeV] 


0.003227 


0.01687 


MSSM 


tan /3 


30 


5 




MA[GeV] 


140 


400 




rA[GeV] 


2.749 


0.7939 




Mh[GeV] 


120.995 


115.948 




rh[GeV] 


0.1047 


0.004083 




MnfGeV] 


140.658 


401.975 




rH[GeV] 


2.640 


0.4726 



Table 1: Higgs-boson masses and decay widths as provided by HDECAY 



the spin-0 channel, i.e. for //^ with equal hehcities, where it appears as a sharp peak at 
s = M^. Apart from the Higgs resonance region, spin-0 contributions to the squared 
lowest-order matrix element |A^oP are suppressed by a factor m^/M| ~ 10~^ with re- 
spect to the spin-1 contributions (/i^ with opposite helicities). This suppression factor 
is clearly visible in the polarized cross sections. At the Higgs resonance an enhancement 
factor M^/Fh in the squared amplitude |A^h,oP compensates this suppression, as long as 
the total Higgs width does not become too large, which is the case if the Higgs decay into 
weak gauge bosons is not yet open. In the latter case (Mh ^ 2Mw), no Higgs resonance 
will be visible anymore over the continuous background induced by the spin-1 contribu- 
tion for realistic degrees of muon beam polarization. Therefore, for the SM we consider 
Mh = 115 GeV and 150 GeV below. 

Now we turn to the photonic and QCD corrections to the polarized cross sections, 
as shown on the l.h.s. of Figure ^ In the spin-1 channels, the cross sections are very 
similar to its counterpart in e"'"e~ annihilation [ |18|, measured at LEP. Above the Z- 



boson resonance, ISR induces very large positive corrections by the so-called radiative 
return to the Z resonance; these corrections even exceed 100% in relative size. The main 
difference between muon and e"'"e~ collisions is due to the different LL factors Asr, which 
are about 0.06 and 0.11 for = M|, respectively, i.e. the ISR effects in muon collisions 
are roughly half as big as in the e"'"e~ case. In the spin-0 channels, which contain the 
Higgs resonance, even much larger positive corrections are visible. These corrections are 
mainly due to a spin flip of the initial-state muons induced by collinear ISR. Note that this 
effect is not accounted for by the LL structure-function approach, which is widely used 
as an approximation for ISR, since this spin-flip correction is not enhanced by a large 
logarithm. The numerical enhancement is rather due to the extremely large difference 
between the spin-0 and spin-1 cross sections. Neglecting (irrelevant) corrections that are 
suppressed by mass factors mj^/s, the spin-flip contribution to the real correction takes 
the universal form \ H7f 



(^-,,spmmp{P+,P-,s) = — / dx(l-x) ao{-P+,P_,xs) + ao{P+,-P_,xs) , (6.4) 

ZTT JO L ^ 

which can be deduced from the factorization of collinear photon emission (see e.g. 



Refs. [^, ^). Thus, the spin-flip correction cr^^spinflip(±l, ±1, s) to the spin-0 chan- 
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Figure 6: Polarized SM cross sections for /i^yU^ bb in lowest order (thin curves) and 
including photonic and QCD corrections (thick curves), where no phase-space cuts are 
applied on the l.h.s. and the invariant-mass cut ( |6.5|) is applied on the r.h.s. 
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Figure 7: Unpolarized SM cross section for [i yU"'' ^ bb in lowest order and including 
photonic and QCD corrections, with and without the invariant-mass cut ( |6.5| ) 
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nel receives large contributions from the lowest-order spin-1 cross sections cro(±l, =Fl) xs), 
multiplied by a/vr and some numerical factor resulting from the convolution over x. Al- 
though these corrections look dramatic for completely polarized beams, they are simply 
part of the large continuous background to the Higgs resonance that is induced by the 
spin-1 channel for realistic degrees of beam polarization.^] 

The large ISR corrections from the radiative return can be suppressed by requiring a 
minimum invariant mass for the bb pair that is sufficiently higher than the Z-boson mass. 
Therefore, we impose the constraint 

y/s - Mhad < 10 GeV (6.5) 

in the following, where Mhad is the invariant mass of all hadrons in the final state, i.e. 
for gluon emission the gluon momentum is fully included in Mhad- This cut implies that 
the energy loss by ISR cannot become larger than 10 GeV. The r.h.s. of Figure ^ shows 
the lowest-order and corrected polarized cross sections in the presence of this cut. The 
corrections to the spin-1 channel reduce to moderate size, as expected, but in the spin-0 
channels still large corrections remain owing to the spin-fiip effect described above. 

The effect of the photonic and QCD corrections on the unpolarized cross section for 
bb is illustrated in Figure |^ for Mh = 115 GeV. Apart from the Higgs resonance 
region, the unpolarized cross section reflects the behaviour of the spin-1 channel discussed 
above in detail. In the following we restrict our discussion to the case of unpolarized muon 
beams {P± = 0), and we always apply the invariant-mass cut (|6.5|) for bb production. 

Figures ^ and || show close-ups of the Higgs resonances for Mh = 115 GeV and Mh = 
150 GeV, respectively. On the l.h.s. the lowest-order predictions for the cross sections are 
successively improved by including the photonic and QCD corrections. The invariant- 
mass cut (|6.5| ) implies that the large positive ISR correction of more than 100% turns 
into a negative correction of ~ —4% (—8%) off resonance and of ~ —49% (—46%) for the 
resonance part for Mh = 115 GeV (150 GeV). This leads to a reduction of the resonance 
peak over the continuous background. The QCD corrections, on the other hand, tend 
to enhance the resonance peak, since they are positive and larger on resonance than 
off resonance. For Mh = 115 GeV (150 GeV) the resonant part of the cross section 
is enhanced by ~ 21% (20%) compared to ~ 5% (4%) in the continuum. This is in 
good agreement with the estimate that is based on the neglect of b-mass effects and of 
interferences between the spin-1 and spin-0 channels. In this approximation the results 



for the decays Z, H — > bb (see Refs. [ |T8|, ^) can be used deduce the following behaviour 
of the QCD corrections, 

A f \ "s(ys) , 17as(yi) 

AaQCD(s) ^ cT^/z,o + ^ t^H,o, (6.6) 

mb^o 71 3 vr 

where ct^/z,o and au,o denote the contributions to the lowest-order cross section ctq induced 
by 7/Z and Higgs-boson exchange, respectively. 

The photonic corrections are dominated by the LL of the ISR corrections in 0{a). 
The remaining photonic contributions are shown on the r.h.s. of the flgures. The largest 

^The relative weight factor between the spin-1 background and the spin-0 signal cross sections is 
roughly given by (1 — P+P-_)/(l + P+P-), which will be significantly larger than a/ir in practice. 
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Figure 8: SM cross section — > bb near the Higgs resonance for Mh = 115 GeV (left) 
and some contributions to the photonic corrections (right): higher-order ISR (h.o. ISR), 
initial-final interferences (IS-FS), FSR, and non-leading 0{oi) ISR corrections (NL 0{a) 
ISR) 
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of these subleading corrections are due to the LL parts beyond 0{a), followed by the 
non- leading (NL), i.e. non-logarithmic, ISR contributions. Note that the latter are the 
part of ISR that is not included in the LL approach via the structure function convolution 



described in Section \i.2.3[ The NL 0{a) corrections reduce the peak over background by 
~ 4% for both Mh values. The FSR and initial-final interaction effects turn out to be 
below 1% of the continuous background and do not show sizable variations in the resonance 
regions. Thus, they can be safely neglected in all phenomenological investigations of the 
Higgs resonance. 

The Higgs resonances inspected in Figures ^ and ^ are so narrow that the beam energy 
spread of the muon collider will eventually determine the resonance shape that is seen 
experimentally. As widely done in the literature, we illustrate this effect by performing a 
Gaussian convolution of the scattering cross section. 



dVI a (VI) exp -^——^ , (6.7) 



where a/s is the central value of the CM energy and the corresponding root-mean- 
square Gaussian spread. Following Ref. [ |T2 



we use 



with the actual value R = 0.003%. The results of applying this convolution to the 
resonances shown in Figures |^ and ^ are shown in Figures |10] and |ll]. As expected, the 
effect of smearing the resonance is much more pronounced for Mh = 115 GeV, where 
the ratio S^/T-n is ~ 0.7, than for Mh = 150 GeV, where the ratio is ~ 0.2. In both 
cases, however, the relative impact of the various radiative corrections remains practically 
unchanged by the smearing. 

6.3 The processes bb, tt in the MSSM 

Now we turn to Higgs production in the MSSM. Of course, it is beyond the scope of this 
paper to give an exhaustive overview over the phenomenology of the full super symmetric 
parameter space. We rather concentrate on the two MSSM parameter choices given in 
Table |l|, which have been selected to be still experimentally allowed and to cover the 
most interesting patterns of Higgs resonances: For tan (3 = 30 and Ma = 140 GeV, the 
lightest Higgs boson h has a mass of M^ ~ 121 GeV with an intermediate width of 
Fh ~ 0.1 GeV, while the resonances of A and H practically lie on top of each other at 
Mh ~ Ma = 140 GeV with large widths Fh ~ Fa ~ 3 GeV. On the other hand, for 
tan f3 = 5 and Ma = 400 GeV, the h boson is SM-like with a mass of Mh ~ 116 GeV and 
a narrow width of Fh ~ 4 MeV; the A and H resonances in this case are close to each other 
but separated, since the mass difference Mh — Ma ~ 2 GeV is sufficiently larger than the 
widths Fh ~ 0.5 GeV and Fa ~ 0.8 GeV. Note also that in the latter scenario the heavy 
Higgs bosons can decay into tt pairs. Since we do not consider the SM-like h boson any 
further, the widths of all resonances discussed in the following are much larger than the 
achievable beam energy spread of a muon collider. Therefore, the smearing hardly affects 
the resonance shapes shown below, and we can omit the convolution ( |6.7| ) in this section. 
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Figure 10: SM cross section fi^ bb near the Higgs resonance for Mh = 115 GeV 
(left) and some contributions to the photonic corrections (right), as in Figure H, but now 
with beam energy smearing 
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Figure 11: Same as in Figure |1^ but for Mh = 150 GeV 
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We first consider the scenario tan /? = 30 and Ma = 140 GeV in Figures 12 and 13 , 
where the h and A/H resonances are shown together with the various radiative corrections, 
in analogy to the previous section for the SM Higgs boson. The photonic corrections are 
dominated by the leading ISR logarithms, and the QCD corrections follow the qualitative 
behaviour of the massless approximation ( |6.6| ), where crH,o is replaced by a"h/H/A,o- Qual- 
itatively the corrections are very similar to the SM expected. The ISR beyond 
0{a) tends to enhance the h and A/H resonances by ~ 7% and ~ 2%, respectively, but 
the NL ISR tends to reduce the resonances by ~ 4%. Photonic FSR corrections stay 
below 1% everywhere, and initial-final interferences do not show a resonance behaviour 
at all. 

Of course, the corrections to the SM-like h-boson resonance of the second scenario 
(tan (3 = 5, Ma = 400 GeV) also do not possess any new feature. Therefore, this resonance 
is not shown explicitly here. More interestingly. Figures |14| and |l^ illustrate the A and H 
resonances, which are close together, for the two processes bb,tt. We note that 

we do not apply any cut like (|6.5| ) in the case of tt production, i.e. the total cross section 
is considered, since a radiative return to the Z resonance via ISR is not possible in this 
channel. For each of the resonances in bb production the structure of the corrections again 
looks as in the previous cases, i.e. the QCD corrections follow the pattern of Eq. 
and the photonic corrections are dominated by the ISR logarithms. The higher-order ISR 
modifies the resonance cross section by ~ 8%, NL ISR by ~ —4%, FSR by less than 
1%, and initial-final interferences not at all. For tt production only minor differences are 
observed. The main differences is the failure of the approximation (|6.6|) , which is due 
to the large mass of the top quark. The leading ISR effects still dominate the photonic 
corrections, while higher-order and NL ISR amount to ~ 6% and ~ —5% on resonance, 
respectively. Owing to the omission of an invariant- mass cut on the tt system, FSR now 
becomes positive, but nevertheless does not exceed ~ 1%, in spite of an enhancement 
by a factor 4 w.r.t. bb production induced by the quark charges. Finally, initial-final 
interferences again do not affect the resonance shape. 



7 Conclusions 



Muon colliders provide a unique framework for studying Higgs-boson resonances, and 
thus for measuring Higgs masses, widths, and various couplings. To this end, it is nec- 
essary to understand the resonance shapes to high precision, i.e. in particular radiative 
corrections have to be controlled. In this paper the photonic and QCD radiative correc- 
tions to — i> // have been discussed in detail, both for the SM and the MSSM. In 
the calculation the issue of a gauge-invariant description of the resonances is emphasized, 
and the corrections that are relevant for the resonance shape are presented in analytical 
form. 

The most important photonic corrections are due to ISR and amount to (9(100%). 
Suppressing the radiative return to the Z-boson resonance by appropriate invariant-mass 
cuts, the calculated ISR corrections, which include the full 0{a) contributions and leading 
logarithmic effects up to 0{a^), should be sufficient to describe photonic corrections within 
per-cent accuracy. Corrections due to photonic FSR modify the resonance by about 1% 
or less. Initial-final interferences (box corrections) turn out to be even smaller and non- 
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Figure 12: MSSM cross section yU^yU^ bb near the h resonance for Ma = 140 GeV and 
tan /3 = 30 (left) and some contributions to the photonic corrections (right) 
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Figure 13: MSSM cross section — > bb near the H and A resonances (which he on 
top of each other) for Ma = 140 GeV and tan /5 = 30 (left) and some contributions to the 
photonic corrections (right) 
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Figure 14: MSSM cross section /i yU+ ^ bb near the H and A resonances for Ma = 
400 GeV and tan /5 = 5 (left) and some contributions to the photonic corrections (right) 
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Figure 15: MSSM cross section yU yU+ ^ tt near the H and A resonances for Ma = 
400 GeV and tan /5 = 5 (left) and some contributions to the photonic corrections (right) 
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resonant, i.e. they modify only the continuous background of the resonance. In summary, 
photonic corrections are now under control at the per-cent level. 

QCD corrections have been calculated to 0{as) including the full mass dependence. 
Further improvements beyond 0{as) are certainly needed in the future. 

The genuine weak corrections are not calculated in this paper. A consistent inclusion 
of these corrections to the full process yU^/i~ —>■ ff (keeping all fermion masses non-zero) 
requires a gauge-invariant Dyson summation of all self-energies, as discussed in some 
detail. However, for a theoretical description of the Higgs resonances only, the relevant 
weak corrections can certainly be incorporated in terms of effective couplings, simplifying 
the task considerably. 
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Appendix 

Standard matrix elements 

In this appendix we give explicit analytical expressions for the SME introduced in 
Eq. (|2.8|) of Section |[ using again the WvdW technique of Ref. [ |^ as in Section |4.2.1 



for the bremsstrahlung matrix elements. Inserting the expressions ( 4.19| ) for the Dirac 



spinors into Eq. (|2.8| ), the SME in terms of WvdW spinor products read 

= {{<P^') + (m') + Mv'Pv)), (A.l) 

where x, y are the vector and axial- vector labels v, a with the sign factors 

= +1, Sa = -1. (A.2) 
Moreover, we have used the shorthands 

1=1,2 i=l,2 



j=l,2 i=l,2 



The explicit expressions for the various helicity configurations are easily obtained by 
identifying the generic spinors 0, ip, etc., with the appropriate insertions according to 
Eq. (pOD. 
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